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We develop a method for the determination of thecdynamics of dissipative quantum systems 
in the limit of large number of quanta N, based on the 1/iV-expansion of Heidmann et al. [ Opt. 
Commun. 54, 189 (1985) ] and the quantum-classical correspondence. Using this method, we find 
analytically the dynamics of nonclassical states generation in the higher-order anharmonic dissipative 
oscillators for an arbitrary temperature of a reservoir. We show that the quantum correction to the 
classical motion increases with time quadratically up to some maximal value, which is dependent 
on the degree of nonlinearity and a damping constant, and then it decreases. Similarities and 
O ■ differences with the corresponding behavior of the quantum corrections to the classical motion in 

' the Hamiltonian chaotic systems are discussed. We also compare our results obtained for some 

, limiting cases with the results obtained by using other semiclassical tools and discuss the conditions 

for validity of our approach. 
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I. INTRODUCTION 



The quantum anharmonic oscillator with the Hamiltonian in the interaction picture (h = 1) 



H = Atfb +J ^(tfb) l+1 , [b,fet] = l (1) 

is one of the simplest and the most popular models describing the quantum statistical properties of light interacting 
with a nonlinear medium 0,^]. In Eq. (|l|), the operators b and b^ describe a single mode of quantum field and the 
constant A; is proportional to the (21 + l)-order nonlinear susceptibility of a nonlinear medium (I is an integer), A 
is the detuning of the light frequency from the characteristic frequency of quantum transition. We adopt the normal 
ordering of operators. For the case of a cubic nonlinearity (1=1), this model first was introduced by Tanas for the 
investigation of self-squeezing of light propagating through a nonlinear Kerr medium without loss. Because of such a 
model is exactly integrable, the explicit time dependence of the quadrature variances necessary for the determination 
of squeezing condition has been found for any moments of time and for any number of photons The problem 

of a dissipative anharmonic oscillator is much more difficult. Nevertheless, Milburn and Holms obtained the exact 
solution for the damped Kerr oscillator (/ = 1) interacting with a reservoir of zero temperature This result has 
been further generalized to the case of a reservoir of non zero temperature in || . 

In the conditions of an experiment, as a rule, a large number of photons arc involved in a nonlinear interaction 
between light and a nonlinear medium modelled by the anharmonic oscillator JljMH] ■ The determination of squeezing 
conditions from the exact solution in this limiting case is straightforward for the model of Kerr oscillator without 
loss In contrast, due to the complex form of the exact solution for the damped Kerr oscillator, the determination 
of photon statistics for the large number of photons in this model demands an application of numerical methods or 
special approximate analytical methods (for a review see 0). Moreover, there are no exact solutions for the model of 
the quantum dissipative oscillator with higher-order nonlinearity and a very little amount of the information on its 
dynamics is documented in the literature |2|]. 

In general, the situation when a large number of photons N are involved in nonlinear interactions is a quite typical 
for many problems of quantum optics |?],|| . Heidmann et al. suggested |J to use the method of 1 /TV-expansion for 
the determination of nonclassical states generation dynamics. They originally applied the 1/iV-expansion technique 
to the problem of squeezing and antibunching of an electromagnetic field interacting with a collection of the Rydberg 
atoms inside a high-Q cavity fiofl , where a large number of atoms is of the same order as number of the photons N. 
The general scheme of the 1 /iV-expansion method states that an exact or an approximate solution of the problem can 
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be found in the classical limit A — > oo and then the quantum corrections could be added pJ| . Because this method 
allows to find the motion equations for the mean values and the lower-order cumulants, it could also be considered as 
a variant of the cumulant expansion [p|JT^] . Recently we further developed and applied the 1 /TV-expansion technique 
Pj to the investigation of an enhanced squeezing at the transition to quantum chaos[] [fu^-|l5|]. It should be noticed 
that only nondissipative quantum systems have been considered in papers |p|,^3|-p^[ . 

In this paper, using the method of 1/A-expansion, we consider a dynamics of squeezing and a deviation from the 
Poissonian statistics for the damped anharmonic oscillators with arbitrary degree of nonlinearity I in the limit of a 
large number of photons N ^> 1. We find the explicit time dependencies for the squeezing and the Fano factor for 
an arbitrary degree of nonlinearity and for an arbitrary temperature of a reservoir. Our consideration is based on the 
quantum-classical correspondence and the fact that the solution of classical equations of motion, obtained within the 
zero-order approximation of X/N, could be found analytically for the case of any linear damping. We show that for 
a weak damping the degree of squeezing is mainly determined by the nonlinear polarization of a nonlinear medium, 
modelled by the nonlinear oscillator. For the case of no damping, our time dependencies for squeezing are transformed 
to the corresponding formulas of work p|, which have been found from an exact solution of the Hamiltonian problem. 
A finite damping decreases the degree of squeezing. The consideration of the Fano factor demonstrates that the 
quantum statistics is always a super-Poissonian for dissipative oscillators. Another restrictive factor having influence 
on the time dependencies of squeezing and the Fano factor are the thermal fluctuations of the reservoir. 

Note that in spite of the fact that we find our main results for the the model of the higher-order oscillator, we 
present our self-consistent system of motion equations for the first- and the second-order cumulants in a form which 
is valid for the description of any single-mode quantum system in the semiclassical limit. One of the main finding 
of this general consideration consists in the influence of a specific quantum diffusive term on the dynamics of the 
expectation values and dispersions. We interpret this diffusion, which is proportional to the damping constant, as an 
influence of the zero-point energy of reservoir on the quantum system. Although the influence of quantum diffusion 
around the classical solution is insufficient for the description of time dependencies of squeezing and mean values 
for the particular system under study, and especially in the most interesting case of a short time of interactions, we 
think that the account of this quantum diffusion is important for the correct description of other dissipative quantum 
systems in the semiclassical limit. 

We compare our basic equations of motion for dissipative systems with the equations arising within the so-called 
generalized Gaussian approximation [I 12 13] and find a one-to-one correspondence up to terms of 1/N 2 for several 



popular models of quantum optics (R - 2l| 

We also discuss the conditions for validity of the cumulant expansion in the form of the 1/A-expansion for the 
description of the dissipative dynamics of nonlinear oscillators. This problem is related to the problem of finding the 
time interval for the quantum-classical correspondence, which attracts large attention nowadays, and especially in 
connection with the studies of quantum chaotic systems (see [ p2| and the references cited therein). For the Hamiltonian 
systems with regular dynamics, the quantum corrections to the corresponding classical equations grow in the course 
of time power- wise ||jj2^,^,^3 15 . As a result, the time interval for the classical description has a power- wise 



dependence on the semiclassical parameter A |2Jj22j,|23[ . In contrast, for the case of nondissipative quantum systems 
which are chaotic in the classical limit, the quantum corrections grow exponentially in time due to underlying local 
instability in the classical system [^4|,^2],^3],[l3] [r| . Therefore, the time interval for the validity of the 1 /A-expansion 
method and the classical description is logarithmic in the semiclassical parameter A |24|,^2|,|3| . 

Our finding for the dissipative nonlinear oscillators is that the quantum correction first increases, then reaches 
some maximum and finally decreases in the course of time. As a consequence, the time interval for validity of the 
1/A-expansion could be divided into two subintervals with completely different dependencies on the semiclassical 
parameter A. The first subinterval scales up power-wise with N and the second one has a logarithmic dependence 
on N. While a power dependence on N originates from the same time behavior as in the Hamiltonian systems 
with regular dynamics, the log A scale has a different nature. It appears because of an exponential damping of the 
underlying classical dynamics. 

Our work is organized as follows. We consider a single mode dissipative quantum system and present the derivation 
of motion equations for the first and second-order cumulants in section [o]. In the same section we also compare 
our approach based on the 1/ A-expansion method with other semiclassical methods and find analytically the time 



1 Recently a good agreement between the predictions for dynamics of squeezing obtained using the method of 1/iV-expansion 
and the results of numerical simulation for the model of kicked quantum rotator with N ~ 10 5 levels were found in the works 
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dependencies for the mean value and the second-order cumulants for the model of the higher-order dissipative oscillator. 
Using these results, we focus on the time dependencies of the squeezing and the Fano factor in section HI. The section 
IV is devoted to the discussion of conditions for the validity of the 1/iV-expansion in the description of nonlinear 
oscillators dynamics. In the concluding section, we briefly summarize our results and outlook the main directions 
for the future developments. Some details related to the solution of motion equations for the cumulants and the 
determination of time interval for the validity of our approach are presented in two Appendices. 



II. BASIC EQUATIONS 

A. The 1 /TV-expansion method and comparison with other approaches 

First of all, we need to generalize the approach of |],0,[l5| to the case of a single-mode quantum system with 
dissipation. Consider an oscillator with the Hamiltonian (|l|) which interacts with an infinite linear reservoir of finite 
temperature. The Hamiltonians of a reservoir and its interaction with oscillator are the following 

H r = ^M d J + V2), H int = (Kjdjtf + H - c . (2) 
j i 

where the Bose operator dj ([dj, d k ] — Sjk) describes an infinite reservoir with characteristic frequencies i[>ji an d Kj 
are coupling constants between the reservoir modes and the oscillator. Introduce new scaled operators a = b/N 1 / 2 , 
cj = dj/N 1 / 2 , and Hermitian conjugats to them with commutation relations 

[a,<f] = l/N, [cj,<*]=6 jk /N. (3) 

In the classical limit N — > oo, we have commuting classical c- numbers instead of operators. Now the full Hamiltonian 
H = Hq + H r + Hi n t may be rewritten as H = NTL, where TL has the same form as in the formulas (jl|) and (||) with 
an account of the following replacements 

b^> a, 6+ — s- a + , dj^Cj, d\ -»• c], and A, -> g t (N) = XN l . (4) 

It could be shown that dependent on photon number constant gi (N) correctly gives the time scale of energy oscillations 
for the nonlinear oscillator (Q) in the classical limit (for the case of the Kerr nonlinearity, see, e.g. 25|). 

Within the standard Heisenberg-Langevin approach the equation of motion has the form ( 0, chap. 7; [ffil ) 



(&-iT)a + V + L(t), (5) 



where V — dTLa/da* , 7 = 2-k\k{<jj)\ 2 p(u}) is a damping constant, p(w) being the density function of reservoir oscillators, 
which spectrum is considered to be flat. The Langevin force operator L(t) — Kjdj(0) exp(— itpjt) has properties 
""261, which in our notations (^) may be rewritten as 

(L(t)) R = (L\t)) R = 0, (it a ) fi+ ( a t J L) fi = 7 M ) ( La ) R + (aL) R = 0, (6) 

where the average is performed over the reservoir variables and {rid} is a mean number of reservoir quanta (phonons), 
related to the temperature T as (rid) = ( c ^(0)c(0)) = [l — exp (t^t)] , where k is the Boltzmann's constant. From 
the Heisenberg-Langevin equations for a, a 2 and Hermitian conjugated equations, we get using Eqs. (|5|) and 

i^-{(6a) 2 ) = 2(V6a) + (W) - iT{(5a) 2 ), (7) 
dr 

ijL/Sa*6a) = -(V*Sa) + (Sa*V) - iT(6a*6a) + iT^$-, 
dr N 

where 

W(a, a*) = ^Tpr, z=(a), ((Sa) 2 ) = (a 2 ) - z 2 , (5a* 5a) = ( a t a) - \z\ 2 , (8) 
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and we have introduced the scaled variables r = git, T = j/gi, A = A/gi. Averaging in Eqs. (Q) is performed over 
both the reservoir variables and the coherent state \a) = exp(Naa' f — Na*a)\0) corresponding to the mean photon 
number ~ N. Such kind of minimum-uncertainty states are most suitable for a consideration of the semiclassical 
limit N ^> 1 pL In derivation of Eqs. (0) we neglect an insufficient additional detuning introduced to A by the 
interaction with the reservoir ELBt|. 

The set of equations (0) is not closed and actually is equivalent to the infinite hierarchic dynamical system for 
the cumulants of different order. To truncate it up to the cumulants of the second order, we make the substitution 
a — > z + 5a, where at least initially the mean z ~ 1 and the quantum correction |<5a(0)| ~ N~ x / 2 <C 1. Using the 
Taylor expansion of the functions V(a,a*) and W(a, a*) around the mean values z and z*, 

— (£),*"(£).*•♦■■■■ —(f)/-©/— o 

and after some algebra analogous to that used in |l^,|f5|, we get from (Q) in the first order of 1/N the following 
self-consistent system of equations for the mean value and the second-order cumulants 

l Tr Z = ~ % \ Z + {V)z + 77° ( Z ' Z *' ((Aa)2> ' ^ Aa *)^' ^ Aa ' 2 ) ' ( 10& ) 
4((Aa) 2 ) = 2^ ((Aa) 2 ) + 2 f|^ N ) (|Aa| 2 ) + {w) z - *r((Aa) 2 }, (10b) 



dr \ da I , V da 



d /IA |2 , fdV*\ //A _, 2x , ( 9V 



i-{\Aa\*) = -^-^-) {{ Aay)+^ — ) ((Aa*)^ - *T(|Aa|") + iT(n d ), (10c) 

where we have introduced the scaled second-order cumulants as ((Aa) 2 ) = N((5a) 2 ) and (|Aa| 2 ) = N(Sa*Sa), as 
well as the function w(z,z*) = NW(z,z*) (for definition of W see Eq. (||)). Now both the first-order cumulant z 
and the second-order cumulants ((Aa) ), (|Aa| 2 ) are of the order of unity a nd sm all parameter 1/N arises only as 
a prefactor for the quantum correction Q to the classical motion equation in ( |l0a[) . The expression for the quantum 
correction has the form of second order differential Q — \d 2 V \ z . In this formula and in Eqs. (|f0|), the subscript z 
means that the values of V and its derivatives are calculated at the mean value z. 



The equations (|l0q), (10c) are nonlinear due to the presence of the nonlinear term (w) z in Eq. (p_0b|). However 



introducing new variables for the second-order cumulants 

B = (|Aa| 2 ) + i CH(Aa) 2 ), (11) 
the term (w) z could be removed and the self-consistent equations (|l0|) may be rewritten as 

ii= (V}, + ±Q(z,z;C,C',B), (12a) 

i6 = 2 (^), c+2 0, s - irc ' (i2b) 
«--.^^-K^). c+ K^)^ + (^).( B -5)- 

If initially the oscillator is in the coherent state, then the initial conditions for system (jlj) are 

5(0) = 1/2, C(0) = 0, (14) 
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and some arbitrary z(0) = zq which is of the order of unity. Involved into equation (12c) an equilibrium value of 



cumulant B is determined by the mean number of reservoir's quanta and its zero-point energy as 

B (0) = M + 1/2- (15) 

Note that the zero-point energy of a reservoir appears in the equations of motion for the cumulants though it was not 
presented in the Hcisenberg equations of motion and even may be dropped from the Hamiltonian redefining a zero 
of the energy. Such "reappearence" of a zero-point field energy is rather often in other problems of quantum theory 
where a vacuum is responsible for the physical effects |27fl . 



It is sufficient that now our basic equations for the second-order cumulants (|12b|) , (12c) arc linear and this fact 
allows us to find analytically the solution of whole set ( |l2|) using the quantum-classical correspondence in the limit 
N —> oo. The substitution ( fil"|) describes, of course, nothing but the transformation from a normal to a symmetric 
ordering. The finding that equations of motion for the lower-order cumulants in the semiclassical limit are solvable 
just within the symmetric ordering is one of the illustration of the fact that symmetric ordering is most suitable for 
the consideration of the quantum-classical correspondence Q . 

Before we will find a solution of equations ([l2]) , let us make several comments and make a comparison with other 
approaches. 

(i) It should be noticed that we derive the equations (^), ( |l2"| ) within the Heisenberg-Langevin approach, but alterna- 
tively, the same equations could be obtained J28| starting from the corresponding generalized Fokker-Planck equation 

§[§■ 

(ii) In the case of no damping T = 0, our equations for the mean values and the second order cumulants ( |l2| ) are 
reduced to the corresponding equations of ref. |[4|,|l5| ■ On other hand, for the Hamiltonians which can be presented 
as sum of a kinetic and a p otential energy, our approach gives the same motion equations as the cumulant expansion 
method introduced in [ p3| . Moreover, Sundaram and Milonni showed p3f that the first-order cumulant approxima- 
tion to the Heisenberg equations of motion gives the expectation values identical to those obtained by the methods of 
Gaussian wavepacket dynamics [ ^9[]3C|| , semiquantum |pl] , ^2[ , and based on the time-dependent variational principle 
|fj3f . In this respect, our present work could be considered as some variant for the generalization of the semiquantum 
methods to the case of a dissipative dynamics. 

(iii) Another approach to the description of the dynamics of quantum fluctuations based on almost Gaussian wavepack- 
ets is the so-called generalized Gaussian approximation (GGA) jl8|,|l],[l2| . This approximation is valid for both the 
Hamiltonian and the dissipative systems and consists in an assumption that the Fourier transform of quantum dis- 
tribution function, i.e. quantum characteristic function, is Gaussian for any moment of time fl2]| . Within GGA only 
the first- and the second-order cumulants are non zero, and therefore it presents the higher-order cumulants in terms 
of only the first and the second-order cumulants and to truncate an infinite dynamic system for the cumulants. We 
compare our system (|l2|) with the corresponding dynamic equations for the mean values and the cumulants obtained 
within GGA for a several popular dissipative models of quantum optics: a second harmonic generation p9| , a non- 
degenerate optical three- wave mixing^] and a forced nonlinear oscillator with cubic nonlinearity pl[ . For the forced 
nonlinear oscillator, we found that our self-consistent set of equations ( |l2| ) coincides with the corresponding basic 
equations of |2l]] up to the terms of the order of 1/N 2 , and for the problems of a nondegenerate and a degenerate op- 
tical 3-wave mixing, our approach gives equations which are identical to the corresponding motion equations obtained 
within GGA. 

There are yet the principal differences between our approach and GGA: First, our approach works good in the 
semiclassical limit N 3> 1, in contrast, GGA is suitable also for the quantum limit N ~ 1. Second, in the general case, 
the equations of motion obtained within GGA are nonlinear and cannot be solved analytically. As we shall see in the 
next subsection, the solution of our system (n2) could be obtained analytically and has simple physical meaning. 



B. The quantum-classical correspondence and the dynamics of cumulants 



We need to find the time dependencies of z, C, and B from the equations (|l2|). First, suppose that we know the 
solution z(t) of the motion equation (12a) for the mean value. Then, the equations of motion for the second-order 



2 The equations of motion within GGA for the model of nondegenerate three-wave mixing could be easily obtained in an 
analogy with the case of second harmonic generation |L9|]. The study of lossless three- wave mixing within GGA has been 
presented in M. 
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cumulants (12b) and (12c) form the linear inhomogeneous equation for the vector variable X(r) = [C(t), C*(t), B(t) 
as 



-i(r/2)X + iX + eX , 



(16) 



where X = (0,0, iTB^) and B^ is given by Eq. (|l5|). The matrix A is formed by the partial derivatives of V 
calculated at the mean value z, its form can be easily obtained from Eqs. ( |l2b| ) and ( |l2c[ ) (to save space we do not 
present an explicit form of A here). We also introduced a parameter e for the discussion of solutions of Eq. Jl6| ) 
(actual value e = 1 in (Eq)). The solution X(r) of linear Eq. (jig) consists of two parts 



X(t) = X(t)+X(t), 



(17) 



where X(r) is a general solution of the homogeneous equation (i.e., Eq. ( |l6[ ) with e = 0) and X(r) is a particular 
solution of the inhomogeneous equation. It easy to see that 



X(r) = [0,TB^T 



(18) 



To find X(r) we need to solve the self-consistent set of equations ( |l 6|) for e = together with (12a) and (l_3|). We 
seek for a solution by the perturbation theory using the smallness of the parameter 1/N. Substituting the expression 



N 



\z (1 \r)\«\z cl (r)\ 



(19) 



into the motion equation (12a) for the mean value, we get in the zero-order of 1/N the classical motion equation 



-i-^Zd + V{z ch z* el ). 



(20) 



Now, it could be shown that solution X(t) of equation of motion for the cumulants (16) with e = can be obtained 
directly from the solution of the classical equation (|2^) by the linearization near z c i (the substitution z c \ — > z c i + 5z, 
\6z\ <§; \z\), if one writes the dynamical equations for the variables (8z) 2 and |<5z| 2 (for the details of derivation, see 
Appendix |X]). Thus, the solution for X(r) is 



X(r)= [(dz cl ) 2 ,(dz* cl ) 2 ,\dz cl \ 2 ], 

where dz c i (r) is the differential of the classical variable z c i (r) governed by Eq. 
are (pT|), i.e. one should take 

(dz ci ) 2 (0) = 0, |d^ c/ | 2 (0) = l/2. 



(21) 

and an initial conditions for (|2l|) 



(22) 



Now we can demonstrate how to find the dynamics of cumulants C and B, if the classical dynamics z c i(t) is known. 
Combining formulas @, @, (HJ) and (§!§), we get 

C=(dz ci ) 2 , B= |dz c/ | 2 + ((n d ) + l/2)rr. (23) 

Turn to the determination of the influence of the quantum correction (|l3|) on the dynamics of the mean value z. 



Substituting the expansion (n9h into equation (12a), we have in the first order of 1/N the motion equation for z^ as 



= + Q(r), Q(r) = Q[^(r), 4(r), C(r), C*(r), B(r)], (24) 

where the cumulants C(r), C*(r), and B(r) are determined by the formulas ( p3| ) and z c ;(r) is a solution of the 
classical equation (pot) . A formal solution of Eq (p3) is 



!l) (r) =z (1) (0)exp^-^r^ -iJ\t'Q(t' 



(25) 



For the initial conditions (14), the initial value of quantum correction ( [13] ) is zero. Therefore, z(0) — z c /(0) = zq and 
z^(0) = 0. The condition of smallness of z^ x \t) in comparison with z c i(t) [Eq. (^9|)] takes the form 
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R(r) 



z{t) - z c i(t) 


1 


z«(r) 


1 


J T dr'Q(r') 


Zci(t) 


~~ N 


Zcl{r) 


~~ ~N 





< 1. 



(26) 



We shall see when the condition (|2g) is valid for the model of anharmonic oscillator (|lj) in the next section. 

We turn to the determination of the dynamics of the cumulants and the mean value in the case of anharmonic 
oscillator with the Hamiltonian ([!]). Here, the expression for V(z, z*) takes the form 

V(z,z*) = Az+ \z\ 2l z. (27) 

The exact solution of the classical motion equation (^o|) with V of the form (^) is 

z(t) = 2 exp [(-iA - r/2)r] exp [-i\z \ 21 ^(t)] , W (r) = [1 - exp(-nr)] /VI. (28) 

Using the expression (Bq), we have from Eqs. (123) the following time dependencies for the cumulants 



C(t) = -m{T)lz%\z \ 2 <- l - 1 '> (»i(r)l\z \ 21 + i) exp [(-r - i2A)r - i2\z \ 21 : w (r)] , 

B(r) = exp(-rr) [l/2 + l 2 \z ftf (r)] + «n d ) + 1/2) Tr, (29) 

where we took into an account the initial conditions (^2|). We shall use these time dependencies for the cumulants in 
the consideration of the squeezing and the Fano factor in the next section. 



III. THE SQUEEZING AND THE FANO FACTOR 



A. The dynamics of squeezing 



Define the general field quadrature as Xg = aexp(—i9) + exp(i9), where 9 is a local oscillator phase. A state 
is said to be squeezed if there exists some 9 such that the variance of Xg is smaller than the variance for a coherent 
state or the vacuum Minimizing the variance of Xg over 9, we get the condition for so-called principal squeezing 




5 = 1 + 2N{(\Sa\ 2 ) - \({Sa) 2 )\) = 2(B - |C|) < 1. (30) 

The determination of the principal squeezing S is very useful because it gives the maximal squeezing measurable by 
the homodyne detection Substituting the expressions (|2^) for the cumulants B and C into the definition (|3C|), 

we find 

S(t) = exp(-rY) [1 + fa(xo, t)} + «n d ) + 1/2)2IY, (31) 



Mxo,r) = 2a a - (1 + a 2 ) 1/2 , a(r) = Ix% 1 ih(t), 



(32) 



where (1i(t) is defined in (|2q ) and we assumed for the sake of simplicity that the initial condition zq is real xq = Rezo- 
First, consider some interesting particular cases. 



1. The weak dissipation 

In the case of weak dissipation Tt <C 1, we have exp(— IY) »1- Tt and Hi(t) « r. Thus, in this limit, from Eqs. 
© and (||), we get 

S(T) = l + (l-rT)<Pi(xo,r) + 2(n d )rT, a(r) « Ix 2 1 t = lT [2l+1) t. (33) 

In the case of no loss (r = 0), the formula ( |33| ) shows that the rate of squeezing is determined by the factor 
21xq L \iN 1 = 2lT'( 2l+1 \ Because of A; is proportional to the (21 + l)-order nonlinear susceptibility, the factor 7- > ^ 2i+1 - ) 
has a physical meaning of the nonlinear polarization. Therefore, the stronger is nonlinear polarization induced by a 
light in the medium, the more effective squeezing of light is possible. For a finite dissipation r ^ 0, the squeezing is 
determined by an interplay between dissipation, the polarization of the nonlinear medium modelled by the anharmonic 
oscillator and the thermal fluctuations of a reservoir. Both dissipation and a finite temperature of reservoir are the 
destructive factors for the squeezing. It is interesting to note that the vacuum term does not appear in Eq. (|3^), 
what is characteristic of the weak dissipation limit. 
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2. The short-time approximation r <C 1 



The short-time approximation r <C 1, as well as the limit of large photon number N 3> 1, are quite realistic for 
a nonlinear medium modelled by the anharmonic oscillators (for numerical estimates, see §, chap. 10, and {§). In 
the limits of r <C 1 and Tt <C 1, we have from Eq. (H|): 4>i(x ,t) ps -Hx^t = -2lV^ 2l+1 h. Thus, substituting this 
expression into (|33|), we get a very simple dependence of 5 on time as 

S(t) = 1 - [lx$ - T(n d )] 2r + 0(t 2 ). (34) 

As follows from Eq. (^), there exists the critical number of phonons n d = (I / ^)V^ 2l+1 ^ such that at rid > an Y 
degree of squeezing is impossible. 



3. The lossless case V = 



In the case of no loss, the time dependence of S is 

S(t) = 1 + 4>i{x ,t), (35) 

where <pi is defined in ( ^2|) and we should take into account that now /i;(t) = r and thus a(r) = Z^q't = lV^ 2l+1 h. 
For the case of Kerr nonlincarity (/ = 1), the formula (|35| ) coincides with the corresponding formula from j| obtained 
in the limit N 3> 1 from the exact solution for the Hamiltonian case. 

Now we consider the time dependence of S [Eq. (|3^)] in the limit of large time r>l, which corresponds to a > 1. 
Rewriting <pi in the form <j>i(xo, r) = 2a 2 [l — (1 + a -2 ) 1 / 2 ] and expanding (1 + a -2 ) 1 / 2 up to the term of order of 
a -3 , we get <j>i k, — 1 + (6a) _1 for a ^ 1. Therefore 

S(t) « (Z^r)- 1 = (V 2 ^ 1 ^ 1 (36) 

for r»l. This is in good agreement with the statement that the principal squeezing S is a power- wise decreasing 
function of time for the general class of integrable systems in the semiclassical limit [|l5|. 

Return to the case of an arbitrary dissipation. The dependence of S on r [Eqs. (|3l|), (|32])1 is shown in Fig. 1, 
where the solid curve corresponds to the lossless case (r = 0), the dashed curve corresponds to the finite damping 
constant T — 5 x 10~ 2 with a reservoir of zero temperature (rid) = 0, and, finally, the dotted curve corresponds to 
the reservoir with finite temperature ((rid) = 1)- As is evident from this figure, the squeezing is stronger for higher 
nonlinearity (compare Fig. la and Fig. lb), a dissipation slows down the rate of squeezing (compare the solid and the 
dashed curves in Fig. 1). Moreover, a finite temperature of a reservoir fastly destroys the squeezing. 



B. The Fano factor 



Another important characteristic of the nonclassical properties of a light is the Fano factor 



F-i^.. ,37) 

(n) 

In Eq. 



which determines the deviation of a probability distribution from the Poissonian distribution with F = 1 
( |37| ) the mean foton number (n) and the mean of square of foton number (n 2 ) are 

(n) = (b%) = N(a?a), (n 2 ) = N 2 (a ! aa t a) = N 2 (a) 2 a 2 ) + (n). 

Substituting expression a —> z + 8a (\5a\ -C \z\ ~ 1) into formulas (|38|), we have after Taylor expansions in the first 
order of 1/N 



(38) 



(a f a) = \z\ 2 + N-\B - 1/2), (a t2 a 2 ) = \z\ 4 + N- l (z*C + c.c.) + N- x A\z\ 2 (B - 1/2), 



(39) 
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and, therefore, the dependence of the Fano factor on the cumulants and the mean values in the first order of l/N is^| 



z 



F = 2B + y—C + c.c.j . (40) 

Substituting expressions ( p9f ) for B and C into Eq. (|l)]), we find the following time dependence for the Fano factor 

F(t) = exp(-rr) + ((n d ) + 1/2) 2Tt. (41) 

As it is evident from Eq. (f4l|), the statistics is the super-Poissonian F(t) > 1 at any time. The most simple form the 
time dependence for the Fano factor takes in the case of weak dissipation 

F(r) = 1 + 2(n d )TT, Ft < 1. (42) 

Thus, the statistics is super-Poissonian for any T ^ and is independent on the degree of nonlinearity This is in 
good agreement with the previous result of ^|J|] for the case of the dissipative Kerr oscillator (I — 1), where from the 
exact solution the impossibility of sub-Poissonian statistics and the antibunching has been predicted. 

IV. THE CONDITIONS OF VALIDITY OF THE 1 /TV-EXPANSION AND THE QUANTUM-CLASSICAL 

CORRESPONDENCE 

The procedure of the 1/iV-expansion may be considered self-consistently, if the influence of the quantum correction 
to classical motion on the dynamics of the mean values is small, i. e. the condition (|2^) is satisfied. From Eqs. jl^ ) 
and (p^), we have the following expression for the quantum correction Q in the case of an arbitrary nonlinearity I 

Q(z, z*) = h(l + l)z* l z l ~ y C + h{l - l)z*V~ 2 h l+1 C* + 1(1 + \)z< l -^z l {B - 1/2). (43) 
We start the consideration of the condition of validity (|2^) with some particular cases. 

1. The short-time approximation 

First, we consider the short-time approximation. In this case, the integral over Q(t) in (^) could be replaced by 
the product and we have 

<1. (44) 

It is easy to see that for r <C 1, Q{t) ~ Q(Q) and thus \(z — z c j)/z c i\ is of the order of l/N. Therefore, the condition 
of validity within the short-time approximation ( [44| ) is fulfilled for any number of photons N 3> 1 large enough. 

2. The lossless case V = 

In Appendix [b| we show that the time interval of validity of our approach for the lossless case T = is 

r«ri =JV i/ 2 |A + M 2 T /2 (45) 

' ^ 'ham JV p/2| ZQ |3/-l - V °> 

This result is in good agreement with our previous finding |l5| that for the Hamiltonian integrable models the time 
scale for the validity of the l/7V-expansion has a power- wise dependence in N in the semiclassical limit N 1. Now 
we turn to the general case T ^ 0. 



z(t) - z c i(t) 


T 


Q(r) 


Zci(t) 


= N 


Zcl(r) 



3 Note that this expression coincides with the corresponding result obtained within GGA (see formulas (3.153) and (10.42) in 
up to the terms of order of l/N 2 . 
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3. The case of an arbitrary dissipation 



The condition of validity fl26|) cannot be obtained analytically for an arbitrary time and for an arbitrary damping T, 
but combining the computations and simple analytic estimates for some limiting case, we can qualitatively understand 
the behavior of R(t) [Eq. (|2(])]. We start with the time dependence of the quantum correction Q(t) [Eqs. (pjj), Jis|) , 
(p8|), (p9|)1- This function is shown in Fig. 2 for the different degrees of nonlinearity (the solid line for I = 1 and the 
dashed line for I = 3), as well as for a rather week (Fig. 2a) and for a relatively strong (Fig. 2b) damping constants. 
The quantum correction Q(t) first increases, then reaches some maximum and finally decreases. The maximum of Q 
increases with the increase of nonlinearity degree (compare the solid and the dashed curves in Fig. 2) and shifts to a 
shorter time with the increase of dissipation (compare Fig. 2a and Fig. 2b). 

Such behavior is in a sharp contrast to the lossless (Hamiltonian) nonlinear systems, where the quantum correction 
always increases with time either power- wise for regular dynamics [ p2[ or exponentially for chaotic dynamics [ p2[ . 
Before Q(t) gets to a maximum, during some time interval [0, t{\ it could be well approximated by the time dependence 



for the lossless case Q{t < t\) ~ 1 3 t 2 (see formula (B8)). If we estimate the time scale t\ from the behavior of the 



most fastly decreasing function ^z(t) [Eq. (pq)], we have 

n = (Tiy 1 , Q( n ) ~ 1 3 t 2 = l/Y 2 (46) 

for the position and the value of maximum of Q(t). As is easy to see from Fig. 2, this naive estimate, however, well 
represents all main features of the time behavior of Q(t). 

Consider now the time dependence of the difference between the classical and the quantum mean values z^'(r) = 
N\z(t) — z c i(t)\ caused by the existence of a quantum correction to the classical motion. Computing the integral over 
Q(t) [Eq. (Pq)1, we plot the time dependence z^'(t) in Fig. 3 for weak (Fig. 3a) and for strong (Fig. 3b) damping, 
as well as for the different degrees of nonlinearity. This plot shows that the difference z^\t) first increases and then 
saturates at some level determined by the degree of nonlinearity and dissipation. 

We start our analysis with the case of weak dissipation (Fig. 3a). For weak dissipation, the level of saturation of 
z^'{t) increases with consequent consideration of higher nonlinearities. This result may be qualitatively understood 
as follows. For weak dissipation, the time interval of almost lossless behavior r < t\ ( f46| ) is rather long and therefore 
the asymptotic behavior of zW (t) sufficiently depends on its behavior during the time < r ^$ T\. For r ^ ri, we 
can use the approximation (BE), i.e. ^^'(ri)) ~ t\I z — l/T 2 . This estimate shows that the level of saturation for 



zW(t) should increase with the growth of I. 

Turn to the case of strong dissipation (Fig. 3b). For strong dissipation, the level of saturation of z^'(t) is sufficiently 
lower in comparison with the case of weak dissipation (compare Fig. 3a and Fig. 3b). Moreover, the level of saturation 
for I = 1 is slightly higher than for nonlinearities with I > 1 and the difference of the saturation levels for the different 
nonlinearities with I > 1 is practically invisible. To understand such behavior we note that for strong enough 
dissipation the time interval for almost lossless dynamics r ;$ n is short and the behavior of Q{t) for another time 
interval \t\, oo] is more sufficient (Fig. 2b). For the nonlinearity with 1 = 1, Q{t) has slowly decreasing tail, while it 
decays rapidly for I > 1 (Fig. 2b). As a result, the area under the curve Q{t) for Z = 1 is greater than the areas under 
the curves corresponding to different I with I > 1. Yet, because all Q(t) are fastly decreasing functions at r > t\ for 
all I > 1, the levels of saturation of z^ x \t) corresponding to different I are indistinguishable in the scale of Fig. 3. 

Till now we studied the time dependence of the difference between quantum mean value and classical solution 
\z{t) ~ z c i{t)\. However, the criterion of validity of the 1/A-expansion in the form ( |2"6| ) includes also the z c /(r) 
as R(t) = N _1 \z^ /z c i\- During the time interval of order of t%, when z^(t) is growing, z c i is slowly decreasing 
and oscillating function. But when z^ x \t) saturates, simultaneously the classical solution should be considered as a 
rapidly decreasing function z c i — exp(— Tr/2). Therefore, for the plateaus in Fig. 3, R(t) ~ TV -1 exp(r-r/2) resulting 
in the corresponding time scale for validity of our approach on a plateau as 

t <t* =r- 1 lniV. (47) 

Formally, this result looks like corresponding estimate for the time scale for validity of the semiclassical description 
("breaking time") in the quantum chaotic systems r* = A _1 lniV, where A is the maximal Lyapunov exponent 



|22| . [23 13-150. However, the physical reason for the appearance of a very short, logarithmic breaking time is quite 
different in our case. When a classical oscillator spiraled around an equilibrium state, the amplitude of the oscillations 
becomes very small and the quantum description starts to differ significantly from classical one due to the uncertainty 
principle. That is the reason for the logarithmically short time scale of the semiclassical description of strongly 
damped oscillations. 
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However, the behavior of an oscillator near an equilibrium point is a physically not very interesting process. On the 
opposite, the physically interesting case when the oscillations are not overdamped, we can again apply the estimate 



(B9) during the time interval r ~ t± = (Tl) and get the criterion of validity in the form 

(1) (rr) 



Zcl(Tl) 



< 1. (48) 



AT 2 



For large enough AT, the criterion Q48J ) could be always satisfied and our approach works well. 

We can summarize our findings concerning the quantum-classical correspondence and the conditions for validity 
of the 1/iV-cxpansion as follows. In a lossless case, the quantum correction to the classical motion equation grows 
quadratically with time. For finite dissipation, the correction Q(t) first increases and then decreases with time. 
The difference between the solution for the quantum expectation value and the classical solution \z(t) — z c i(t)\ first 
increases a power-wise with time and then saturates at some level, which is dependent on the degree of nonlinearity 
and on the dissipation strength. The condition for validity has different scaling dependence on the semiclassical 
parameter iV for two different time intervals. During the time interval of growth of \z(t) — z c i(t)\, the criterion of 
validity has the form fl48|), i.e. R ~ X/N <C 1, and it is certainly satisfied for large N. During the time interval when 
\z(t) — z c i(t)\ saturates, our approach is valid only for a short time t* ~ In AT. We believe that described above 
conditions for the semiclassical description of the quantum system are common for all dissipative systems with simple 
attractor. 



V. CONCLUSION 



In summary, using the quantum-classical correspondence, we find analytically the time dependencies of squeezing 
and the Fano factor for the dissipative anharmonic oscillators of an arbitrary degree of nonlinearity in the limit of 
large number of fotons. It should be noticed that our basic equations are rather general and could be considered as 
some generalization of the Gaussian wavepacket dynamics methods to the dissipative systems. Moreover, some effects 
mentioned in the present work are also rather general. For instance, quantum diffusion around classical trajectory in 
the dissipative system, which we interpret as the influence of the reservoir's vacuum on the quantum system, should 
be observable not only in the dissipative nonlinear oscillators. In this respect, we would like to mention the work 
of Savage published ten years ago Q]. The author of |54j studied numerically, within the Gaussian approximation, 
the quantized version of the second harmonic generation problem for the self-oscillating regime 35 1. He had found 
numerically a diffusive growth of the quantum variances near a classical limit cycle. However, neither an analytic 
treatment of the problem, nor a physical explanation of this "diffusion" were presented in J34]. We plan to devoted 
our future publication to the investigation of this problem together with more accurate study of the validity conditions 
for the application of the l/A^-expansion to the dissipative quantum systems. 
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APPENDIX A: 

In this Appendix we show that the equations of motion for the cumulants (|l^) at e = could be obtained from the 
classical motion equation ( pp| ) . This proof is valid for any dissipative system with one degree of freedom. 
Linearization of classical equations ( p0[ ) near z c \ by means of the substitution z c i — » z c i + Sz gives 

d r dV dV d T dV* dV* 

i—Sz — ~i—Sz + ——Sz + — — Sz*, i—5z* = —i—Sz* — Sz — „ — Sz*, (Al) 

dr 2 oz oz* dr 2 oz Oz* 
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where all derivatives are taken on the classical trajectory z c i(t). Using the equality (dV/dz) = (dV* /dz*), we have 
from Eq. (Al) for the quadratic variables (Sz) 2 and |<5z| 2 the following equations of motion 



9 dV 9 dV 
-iT(6zf + 2^(Szf + 2^Sz 2 , 
oz oz* 

T\X 12 dV * (X \ 2 . 9V fX *\ 2 

-iTSz -——(Sz) + — (Sz) . 
oz oz* 



(A2) 



It is easy to see that Eqs. (A2) are equivalent to Eqs. ( |16|) with e = 0, if one makes substitutions \Sz\ 2 — > B and 
(Sz) 2 -» C. 

It should be noticed that yet exists one difference between the linearization of the classical motion equations and 
the equations for quantum cumulants (12b), (12c): It is impossible to get the initial conditions (14) for C and B from 
only initial conditions for the linearized classical equations of motion (see also discussion of this problem in fll4yi5| ) . 



APPENDIX B: 



In this Appendix, we present the derivation of the time scale for the validity of the X/N- method for lossless (r = 0) 
oscillators. Start with the case of Kerr nonlinearity I = 1 and then generalize obtained results for arbitrary I. From 
Eqs. ( p8| ) and (29) for V — and I = 1, we get 

Zd{r) = z exp(-iOr), fl = A+\z Q \ 2 , 



C{t) = -rzl [\z \ 2 t + i] exp(-i2fir), B(t) = 1/2 + |z | 4 r 2 . 
The expression ( |43| ) for quantum correction in this case is 

Q(z, z*, C, C*,B) = z*C + 2zB - z. 



Substituting ( JBl| ) and ( p32[ ) into (p33|), we have for defined in Eq. function Q(t): 

Q(t) = 2;o|zo| 4 ''" 2 exp(— iilr) — izo\zo\ 2 T exp(—iQT), 

and involved in ( p6| ) integral is 

/ dr'Q(r') = z ko| 4 ^ 3 [2fiTexp(-ifir) + i{fl 2 T 2 - 2) exp(-zfir) + 2i] 
Jo 

izo\zo\ 2 Q~ 2 [exp(— iilr) + ifirexp(— iVLr) — 1] . 



(Bl) 
(B2) 

(B3) 

(B4) 



(B5) 



The most rapidly increasing term in (B5) for r 3> 1 is f2 ^^zol^o^r 2 exp(— if2r). Substituting this expression into Eq. 
(Eq) , we have condition of validity in the form 



N 4 r 2 

\A+\z \ 2 \ N 



«1, 



and therefore the time scale r* of validity is 



r<<T . = „v,iA + w'r 



(B6) 



(B7) 



Turn to the lossless case for an arbitrary L Analogously to the case I = 1, it could be shown that the expression for 
quantum correction (|24|) in the limit of large time r 3> 1 is 

Q(r) w zoNol 2 ^" 1 ^ 3 ^ exp[-tn,r], fi z = A + |z | 2i , (B8) 

and the influence of quantum correction on mean value (^q) is 



2 (D(r) = / dr'Q(r') 



^ 3 ^ |zo| 2 ( 3 ^ 1 )exp(-zr!r^ 



(B9) 



where we again presented only the term which is leading in time for t»1. Finally, substituting (|B9j) to (|26|) with 
account of (Eq), we obtain (45). 
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FIG. 2. The dependence of the quantum correction Q on scaled time r for weak damping Y = 0.05 (a), and for strong 
damping V — 0.5 (b). Solid curve corresponds to the Kerr nonlinearity with I = 1 and dashed curve - to nonlinearity with 
I — 3. Everywhere (rid) = 0, A = 0, zo = 1. 
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FIG. 3. The scaled difference between the quantum and the classical expectation values \z^\ = N\z — z ci \ plotted as a 
function of time r for weak damping V = 0.05 (a), and strong damping Y — 0.5 (b). The Kerr nonlinearity with I — 1 is 
described by the solid curve, dashed curve corresponds to the nonlinearity with I — 3, dashed and dotted curve corresponds to 
the nonlinearity with I — 5. Other parameters and initial condition are the same as in Fig. 2. 
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